Phase diagram and effective shape of semi-flexible colloidal rods and biopolymers 
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We study suspensions of semi-flexible colloidal rods and biopolymers using an Onsager-type 
second-virial functional for a segmented-chain model. For suspensions of thin and thick fd virus 
particles we calculate phase diagrams in quantitative agreement with experimental observations, 
and we find their effective state-point dependent shape to be much shorter and thicker than the 
actual shape. We also calculate the stretching of worm-like micelles in a host fd virus solution, again 
finding agreement with experiments. For both systems, our results show that the fd virus stiffness 
can play a key role in system behavior. 



Rod-like particles are capable of forming a great vari- 
ety of liquid-crystalline phases [l| , and have been widely 
explored experimentally and theoretically 0]- Recently 
a prominent role is being played by aqueous suspensions 
of fd virus particles, which are charged, semi-flexible, 
Brownian needles with a length-to-diameter ratio exceed- 
ing 100, exhibiting isotropic (I), cholesteric nematic (N), 
smectic, columnar, and crystalline phases upon increas- 
ing the concentration 0-0]. Moreover, wild- type fd virus 
particles have been bio-engineered to have, for instance, 
a polyethylene-glycol (PEG) coating, such that mixtures 
of thin and thick rods with diameter ratios d varying from 
1.1 to 3.7 could be studied experimentally @. The re- 
sulting experimental phase diagrams are extremely rich, 
even in the regime where only I and N phases are rele- 
vant. For instance, the observations not only include I-N 
coexistence with strong fractionation effects, but also, for 
d ^ 3, two-phase N-N and three-phase I-N-N coexistence 
with phase diagram topologies that strongly depend on 
d. Given the needle-like shape of the fd virus particles, 
and the relative structural simplicity of I and N phases, 
one would expect these phase diagrams to be well un- 
derstood, e.g. resembling those of theoretical predictions 
based on Onsager's second virial theory for thin-thick 
mixtures of rigid rods 0- However, these systems turn 
out to be surprisingly poorly understood. For instance, 
the experiments show N-N demixing with a low-density 
(lower) critical point that shifts to lower densities with 
increasing d @ , while the theory predicts the exact oppo- 
site: a high-density (upper) N-N critical point and N-N 
demixing that extends to higher densities with increasing 
d Moreover, the experiments show N-N demixing at 
diameter ratios as low as ci ^ 2 — 3 while rigid-rod mod- 
els predict d ^ 4 — 5 [1, Q • In this Letter we will show 
that the key to a real understanding of these systems is 
flexibility, which renders the needles effectively shorter 
and fatter depending on the state point [1, 0] . In addi- 
tion, our model and theory can also quantify the observed 
stretching of guest biopolymers in host suspensions of fd 
virus particles [IQ] . 



One-component systems of semi-flexible rods have 
been studied using numerous methods, and it is known 
that only a slight flexibility is enough to shift the I-N 
transition to significantly higher concentrations [lll - [l5| . 
We build on the segmented-chain model introduced by 



Wessels and Mulder [15|, in which (i) flexibility is incor- 



porated by introducing a bending potential between the 
chain segments, and (ii) excluded volume is taken into 
account at the segment level. This approach reproduces 
(in the appropriate limits) the results of Ref. pjj], with 
the advantage of only having to deal with a discrete num- 
ber of degrees of freedom, making explicit calculations of 
extensions to two-component systems feasible. 

We consider a suspension of Ni semi-flexible rods of 
species i — 1,2 with contour lengths , in a volume V 
at temperature T. Following Ref. [l^ we model a rod of 
species i as a chain of Mi rod-like segments of length U = 
Li/ Mi and diameter Di <C U- Denoting the orientation 
of the m-th segment by a unit-vector (with 1 < to < 
Mi), we write the bending energy of a chain of species i 
with orientation $7 = {cji , . . . , wa/. } as 



m— 1 



(1) 

where the stiffness is described in terms of the persistence 
length Pi [15^. Here and below we use thermal energy 
units by setting ksT = 1. The state of the suspension 
is characterized by the orientation distribution functions 
(ODFs) fi{^), which satisfy the normalization condition 
J drtfi{fl) — 1 where dft — Y[m=i'^^m- Denoting the 
total number of rods by TV = N1+N2, the density hy p = 
N/V, and the mole fraction of species ihy Xi = Ni/N, we 
can write the variational free-energy functional F[fi, f2\ 
within an Onsager-like second virial approximation as 
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ln{Bp) — 1 + xi In xi + xi In xi 
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The first line of Eq. ([5]) represents the translational and 
the mixing ideal-gas contributions (with B = ^DiL\, a 
constant), the second line, the orientation entropy and 
bending energy, and the third line the excluded volume 
interactions, which can be considered at the segment level 
with Kij{Vl, Jl') = I]mli SmLi hjiujm, Wm')- The free- 
energy functional in Eq. ([8]) is a two-component general- 
ization of the one-component segmented-chain functional 
of Ref. [HI, and for Mi = 1 and Ui = it reduces to 
the Onsager functional for binary mixtures of rigid rods 
0j[3- From hereon, we simply summarize our method, 
and direct those interested in the complete outline to the 
appendix. At a given thermodynamic state point, the 
equilibrium ODFs minimize F and therefore satisfy 



nm = 



(3) 



v,{n) = pY,xj j K,,{n,n')f,{n')dn\ (4) 

where Vi{fl) can be seen as a self-consistent field act- 
ing on all segments of a chain, and Qi is a partition 
function-like normalization factor Explicitly solving 
Eqs. ([TU)) and ([TT]) for state points of interest would be 
prohibitively expensive computationally because of the 
high-dimensional angular il-grids that would be required 
in the case when Mi ^ 1. Instead, we formally evaluate 
the functional F of Eq. ([8]) in its minimum by inserting 
the solutions fi of Eqs. (fTOl) and ([TT]) to find the equilib- 
rium free energy 

^ = \n{Bp) - 1 + 2^1 In ^ + 0:2 In ^ - ^ x^Xj 
X ^ ^ / kij{uj,uj')fi^jnioj)fj,m'{uj')dujduj', (5) 

m=l va' = l 

where fi^m{i^) is the ODF of the m-th segment (to = 
1, . . . ,Mi) of a chain of species i = 1, 2 given by /i,m(ci;) — 
J fi{S},)5{uj—u!m)dfl. Eq. ([12]) implies that the thermody- 
namics does not require the full solutions fi{^) but only 
the Mi single-segment distributions fi^,n{uj) and the nor- 
malization factors Qi. The calculation of these is given 
in the appendix. With the ODFs, and hence Feq, known, 
all thermodynamic properties such as osmotic pressure 11 
and phase diagrams follow Q- 

We have calculated the phase diagrams for mixtures 
of bare fd particles (species 1, thin) and PEG-coated 
ones (species 2, thick), with equal contour and persis- 
tence lengths, Li = L2 = 0.88/im and Pi — P2 = 2.2/ini 



[6j. The bare fd diameter is fixed to Di = 6.6nm, and 
following the experiments of Ref. Q we consider several 
diameter ratios d = D2/D1 to describe varying thick- 
nesses of the PEG coating. Throughout, we use suffi- 
cient segments per virus particle such that we are in the 
continuum limit for all state points of interest (see the 
appendix). In Fig. [1] we show, for several d, the resulting 
phase diagrams in the (»7i,?72) representation, where rji 
denotes the packing fraction of species i — 1,2. For all 
d we find isotropic-nematic (I-N) coexistence, with the 
tie-lines indicating, for increasing d, an increasing frac- 
tionation of the thinner and thicker rods preferentially 
into the I and N phase, respectively. 
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FIG. 1. Phase diagrams (see text) for mixtures of bare (thin) 
fd virus particles (species 1) and PEG-coated (thick) ones 
(species 2), with diameter ratios (a) d — 1.2, (b) 3.1, (c) 3.11, 
(d) 3.125, (e) 6. The lighter colored areas indicate the two- 
phase regions with tie-lines connecting coexisting state-points; 
triangles denote I-N-N and I-I-N triple points, (f) Topology 
of the phase diagram (two- and three phase coexistence) of 
binary mixtures of (modified) fd virus particles as a function 
of their diameter ratio d and persistence length P. 

In agreement with the experiments of Ref. [6] we find 
no N-N coexistence in the density regime of interest for 
the smallest diameter ratio d = 1.2 (Fig. (Ua)). For in- 
creasing d an N-N demixing regime appears in the density 
regime of interest, with a lower critical point as shown 
for d = 3.1 in Fig. [T](b); in the experiments of Ref. @ 
such a phase diagram was found for d ~ 2.9. Slightly 
increasing the PEG layer thickness to d = 3.11 reveals 
an I-N-N triple point and an associated N-N demixing 
regime with an upper critical point emerging out of the 
I-N coexistence regime (Fig. [T](c)). This N-N upper crit- 
ical point has not been reported experimentally, perhaps 



3 



because it only exists in a small regime of diameter ra- 
tios: for d = 3.125 (Fig. [Ijd)) the upper and lower N-N 
points have merged to form a single neck-shaped N-N 
regime that is separated from the I-N coexistence regime 
by an I-N-N triple point, with the I and one N phase 
rich in bare fd-particles and the other N phase rich in 
PEG-coated fd-particles. Interestingly, such neck-shaped 
phase diagrams have also been reported in Ref. :0] for di- 
ameter ratios d ~ 3. For d — 6, Fig. [Ije) shows an 
I-I-N triple point and I-I coexistence, as found also for 
thin-thick mixtures of rigid rods with diameter ratios ex- 
ceeding d ~ 8 0. We find that this behavior is present 
for d > 4.5. 

Motivated by recent progress in the bio-engineering of 
fd virus particles [l7| . which may allow for tuning their 
flexibility, we have also calculated binary-mixture phase 
diagrams for a large variety of diameter ratios and per- 
sistence lengths Pi = P2 = P. Fig. [TJf) summarizes our 
flndings by dividing the (d, P) plane into regimes with 
phase diagrams featuring only an I-N transition (for stiff 
rods and small d) all the way to complex phase diagrams 
with I-N, I-I, N-N phase coexistence and I-I-N and I-N- 
N triple points (for flexible rods and large d). Clearly, 
increasing d and decreasing P have similar effects on 
the phase diagram, and hence increasing the flexibility 
is expected to considerably enhance the complexity of 
the phase diagrams. 

In order for a model system of rigid rods to fully cap- 
ture the phase behavior of binary fd virus systems, one 
must use a system of much shorter, thicker rods, with 
L/D ^ 7 for the PEG-coated fd virus Q, much lower 
than the true values of 20 — 110. The inference is that 
long semi-flexible rods exhibit the same phase behavior 
as short rigid rods [1, 0] . Our model also enables us to 
study the effective shape exhibited by the flexible rods 
throughout the binary fd virus systems. We define the 
mean square effective length ^ as 

Mi Mi 

^1,1 = X! X! / ^'^ ' ^')fi,m,m'{^,^')dujduj', (6) 
m— 1 m' — 1 

where fi^m,m'i'-^,'-^') is the pair orientation distribution 
function (PDF), the calculation for which is given in the 
appendix. We use the PDFs to calculate Le,i and, from 
this, the effective diameter D^j required for rigid rods of 
length Le,i to have the same excluded volume as flexible 
rods of length Li. 

Fig. [2ja) shows the effective shape of the rods 
Li,s/De,i, for a mixture of thick-thin fd virus particles 
with d = 3, throughout the phase diagram (the route we 
follow is shown in Fig. [UJb) , where we use the X2 — ^ 
representation for clarity). What is immediately ap- 
parent is that throughout the phase diagram, while the 
rods always behave as shorter, thicker rods, the effective 
shape varies considerably. In the isotropic phase, we flnd 
Le^i/De^i to be about 20% of Li/Di (for both species). 
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FIG. 2. (a) Relative effective shape [Le,i/ De,i)* ~ 
{Le,i/ De,i) / {Li/ Di) of fd virus particles as a function of os- 
motic pressure (11) (solid lines show bare fd, dashed lines show 
PEG-coated fd), and (b) the path we follow throughout the 
phase diagram in the X2 — ^ representation, indicated by 
the arrows. 

This corresponds to Lf./Dc. w 8.5 for the thick rods, close 
to the L/D ^ 7 required for rigid rod systems to capture 
the phase behavior of these binary systems 0. For the 
nematic phase, however, L^^i/De^i jumps to over 50% of 
Li/Di, and increases considerably as 11 is increased. We 
conclude that a flxed effective shape does not capture the 
essential physics of these suspensions; the state-point de- 
pendent stretching of the flexible rods is a key feature. 

Finally, we present results for the effective length of 
semi-flexible polymers dissolved in an fd virus suspen- 
sion. A range of polymers which undergo a coil-rod tran- 
sition, stretching out over the I-N transition of the host 
fd virus, has been studied experimentally [lOj. Here, we 
examine worm-like micelles, which have constant P = 
0.5/im, D = 14nm and variable L = 5 — 50/im. The con- 
centration of the polymers is suflSciently low that they 
can be treated as a single particle in a bulk fd virus sus- 
pension, and we study the behavior of the polymer L^. 
over the I-N phase transition of the fd virus (which is at 
UB = 29.54). The results are shown in Fig.|31[a). For the 
shortest polymers studied, we see a considerable jump in 
Le, from Lf. ~ 0.26L to ~ 0.61L, corresponding to a 
coil-rod transition. For longer polymers, the jump in Lf. 
becomes smaller, and the longest ones only become truly 
rod-like well into the nematic phase of the fd virus. It 
is interesting to note that in the isotropic phase, for all 
cases, Lf. appears to remain essentially constant. 

In the isotropic fd virus phase, the polymer may be 
considered along the lines of the Kratky-Porod worm- 
like chain mode l Il8i , generalized to account for excluded 
volume effects [l^- The average end-to-end length is 
deflned as Lkp = V4P'P(L/2P)'', where P' is an effec- 
tive persistence length. For an ideal Kratky-Porod chain 
P' — P and V = 0.5. We calculate L^ in the isotropic 
phase (Lej) for a large range of L values, and flt Lkp to 
our results for the range of worm-like micelle parameters. 
We find P' = 0.573P, and v = 0.529, shown in Fig.EJb), 
where we also compare our results to the ideal Kratky- 
Porod worm- like chain, and to the rigid rod length L. 
Glearly, for shorter micelles, Lf.j approaches L, whilst 
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FIG. 3. (a) Le of worm-like micelles of various L against host 
fd virus osmotic pressure 11. (b) Le.i against L, with fit to 
Lkp, compared to the ideal Kratky-Porod Lkp and rigid rod 
L. (c) Le,N of worm-like micelles with L = 10/^m, and Lo 
of a confined semi-flexible polymer, against 11. (d) AL^ of 
worm- like micelles against L/P of host fd virus, compared to 
AL' 



for longer ones, L^j approaches the ideal Kratky-Porod 
Lkp- 

In the nematic fd virus phase, we can consider the poly- 
mer using the Odijk confined worm-like chain model 14| . 
where the host nematic acts as a confining cylinder. The 
average end-to-end length is defined as Lq — L{uj ■ n), 
proportional to the average lu ■ n along the chain, where 
n is the nematic director. We find that our calculated 
effective length at the TN phase transition L^^n is sig- 
nificantly below Lo (Fig. injc)), from which we infer that 
the host nematic is neither dense nor ordered enough to 
fully confine the polymer. By altering the stiffness of 
the host fd virus such that it is denser (flexible fd virus) 
or more ordered (rigid) at the transition (see the ap- 
pendix), we may confine better the polymer causing it to 
stretch out more over the I-N transition. Fig. ^d) shows 
ALe — Le^N — Lej for worm-like micelles of L = lO/xm, 
using a range of P values for the fd virus. Here we see 
that ALe is the largest when using the most flexible fd 
virus and it decreases as the fd virus becomes more rigid, 
until it reaches a minimum, before increasing again for 
the most rigid fd virus studied. This presents the oppor- 
tunity to tune the stretching of the polymers by varying 
the stiffness of the host fd virus. We also compare ALe 
to Aig ~ Lo — Lkp (Fig. [3Id)), flnding that this over- 
estimates ALe, but does qualitatively match our non- 
monotonic results. 

In conclusion, we have developed a model for binary 
mixtures of semi-flexible rods, and applied it to binary fd 
virus mixtures. We find I-N, N-N (with upper, lower and 
no critical points) and I-N-N coexistence regions present 
in such systems, consistent with experimental results 
Additionally, we have shown that the N-N upper critical 



point only exists in a very small diameter ratio regime, 
which may explain why it has never been seen experi- 
mentally. For systems with diameter ratios larger than 
those studied experimentally so far, we also find TI and 
I-I-N coexistence regions. We also see that altering the 
stiffness of the rods can have the same effect on phase 
behavior as altering the diameter ratio, giving an extra 
parameter to tuning the phase behavior of rod-like parti- 
cles. The addition of flexibility in our model gives results 
that are quantitatively closer to experimental observa- 
tions than those obtained using rigid rods Q. The key 
reason appears to lie in the effective shape of the rods, 
which changes throughout the phase diagrams, indicating 
that any static-shaped rigid rod model will miss some of 
the essential physics. We have also studied the stretching 
of semi-flexible polymers in an fd virus solvent. We find 
that sufficiently short polymers stretch out considerably 
over the I-N transition of the host solvent, while for the 
longer ones, some stretching is observed, but the effect 
is much less pronounced. Changing the stiffness of the 
host fd virus, such that the density or nematic order at 
the I-N transition is increased, can greatly increase the 
stretching effect. We hope that our findings will stimu- 
late further experimental explorations, and believe that 
extensions of the theory to e.g. inhomogeneous states 
and non-equilibrium phenomena are within reach of the 
present model. 
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APPENDIX 
Theory 

Free-energy minimization and phase diagrams 

In the Onsager theory (20| , the isotropic- nematic phase transition of a one-component system of rigid- rods is driven 
by competition between two entropies. One is similar to entropy of mixing, arising from the mixing of particles of 
different orientations. The other arises from the excluded volume interactions of the particles. These are expressed 
via the Helmholtz free energy, F, and for monodisperse rigid particles the solutions to this are well known (see e.g. 



2lLl22|). Extending the theory to polydisperse systems requires the addition of an extra mixing term to the Helmholtz 
free energy 0, [3]- Flexibility may be incorporated in numerous ways (see e.g. 11- IE]), and in this supplementary 
material, we build on the work of Wessels and Mulder [l^ to produce a model describing binary systems of semi-flexible 
rods. 

We consider a suspension of semi- flexible rods of species i = 1,2 with contour lengths Li, in a volume V at 
temperature T. Following Wessels and Mulder [l^ we model a rod of species z as a chain of Mi rod-like segments of 
length li = Li I Mi and diameter Di <^ li. Denoting the orientation of the m-th segment by a unit- vector Um (with 
1 < m < Mi), we write the bending energy of a chain of species i with orientation Q, — {loi, . . . , to Mi} as 

A/i-l p Mi-1 

Ui[Q.) ^ ^ U.i{LOm,UJra+i) = X! ^rn ' (7) 

iTL—l ^ m — 1 

where the stiffness is described in terms of the persistence length Pi ^15i] . Here and below we use thermal energy units 
by setting fc^T = 1. The state of the suspension is characterized by the orientation distributions functions (ODFs) 
/i(r2), which satisfy the normalization condition J dflfi{ft) — 1 where dfl = rim=i dj^m- Denoting the total number 
of rods by iV = A^i -I- 7V2, the density hy p = N/V, and the mole fraction of species ihy Xi = Ni/N, we can write the 
variational free-energy functional F[fi, /2] of this system within an Onsager-like second virial approximation as 

F 

— — ln{Bp) - 1 + xi In + a;2 In X2 (8) 



N 



2 „ 

/ /.(fl)(ln(47r/,(n)) + f/,(n))dn 

J2 ^^^3 J f^{^)f](P■')K^A^,^')d^dn'. 



The first line of Eq. ([5]) represents the translational and the mixing ideal-gas contributions (with B — ^DiL\, 
a constant), the second line, the orientation entropy and bending energy, and the third line the excluded volume 
interactions given by 

Mi Mj 

Kij{p.,0.') = ^ X! kij{uJm,UJm') 
m=l m' = l 

Mi Mj 

= lilj{Di+Dj)^ ^ I sin7(cj,m,cj„/)|, (9) 

m— 1 m'—l 

with 7(wm,ajm') — arccos(a;m • oJm') the angle between chain segments [l5| m and m'. The free-ener gy f unctional 



of Eq. ([5]) is a two-component generalization of the one-component segmented-chain functional of Ref. 15|, and for 
Mi = 1 and Ui = it reduces to the Onsager functional for binary mixtures of rigid rods [tI. [l6|. 

At a given thermodynamic state point, the equilibrium ODFs minimize F and therefore satisfy the Euler-Lagrange 
equations S{F — iJ,iNi)/dfi{^l) = for i = 1,2, with fii the chemical potential-like Lagrange multiplier that ensures a 
proper normalization. This gives rise to 
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= pY,x3 f K,j{n,n')fj{n')dn\ (ii) 

where Vi(r2) can be seen as a self-consistent field acting on all segments of a chain, and Qi is a partition function- like 
normalization factor. Explicitly solving Eqs. (jlOp and (jlll) for state points of interest would be prohibitively expensive 
computationally because of the high-dimensional angular fJ-grids that would be required in the case when Mi ^ 1 . 
Instead, we formally evaluate the functional F of Eq. ([S]) in its minimum by inserting the solutions /,; of Eqs. (jlOp 
and to find the equilibrium free energy 

= ln{Bp) - 1 + xi hi — !- 4- 2:2 In - -p Xia;^ / % (w, uj')fi.m (i^)/j,m' {uj')dujduj' , (12) 

I J m=l m' = l 

where fi,m{^) is the ODE of the m-th segment (m = 1, . . . , M^) of a chain of species i ~ 1,2 defined by 

fi,m{uJm) = y" fi{^)duji . . . duJm^idujjn+1 ■ ■ ■ dcjMi- (13) 

Eq. (IT^ implies that the thermodynamics does not require the full solutions fi{fl) but in fact only the single- 
segment distributions /i,m(w) and the normalization factors Qi, for which an efficient iterative recursion scheme, that 
exploits the connectivity of the chain, can be set up as follows. 

Eqs. ([9]) and ((TT|) allow us to write Vi{Ct) — X]m=i ^'ii'^m) with the same selfconsistent field 

2 Mj 

,{uJm,UJm')f],Tn'{,UJm')duJm', (14) 



,(w™) = Xj / h 

1 = 1 m' = l 



for all segments of chains of the same species. As a consequence, Eq. ([TO)) combined with Eq. (fT3)) can be written as 

/j,m(w) = -^9i,„(w)exp[-Wj(a;)]gi3f_,„+i(iLj), (15) 
with the partial-chain partition function 

/ 771—1 
exp[-Ui(a;„) - Ui(iu„,tJ„+i)](iiu„. (16) 

77=1 

In the formulation of Eq. p3|) the m-th segment ODE is seen as the statistical weight exp(— Wi(aj)) of that segment 
in the (selfconsistent) field Vi{uj), combined with the weights qi^m and qi^M-m+i of the two sub-chains from segment 
m ± 1 to the two chain ends, respectively. Interestingly, the connectivity of the chain allows us to rewrite Eq. (1161) as 
the recursion relation 

97,777(1^) = J qi,m~iiu}')exp[-Vi{uj') - Ui{uj',uj)]du;' (17) 

such that a loop can be set up that (i) starts with a guess for Vi{uj) for i = 1,2, (ii) solves for qi^m{'-^) for all i = 1,2 
and m = 1, . . .Mi using Eq. ([T7|) together with qi^iiui) = 1, (iii) computes fi^m{i^) and the normalization factor Qi 
from Eq. (fTS]). (iv) recalculates Vi{oj) using Eq. (fT4|) and repeats (ii)-(iv) until convergence is found. Note that this 
scheme only requires an angular grid for uj, which in the light of the azimuthal and up-down symmetry of the isotropic 
and nematic phases of interest here, reduces to a single grid for the polar angle £ (0, 7r/2). The only difference 
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with Onsager-type theories for rigid-rod mixtures is the additional calculation and storage of for 1 < m < M, 

here. With the ODFs known, Feq may be calculated from Eq. (IT2|) . The osmotic pressure is then calculated from 

n = P^^~^r~- ^"-"^ ^ binary system, phase behavior is most easily analyzed using the Gibbs energy per particle, 
g{x, n) = + — . By fixing 11, g may be calculated as a function of X2 (with xi = 1 — 2:2), and performing a common 
tangent construction allows for the prediction of coexisting phases 0, 0, [lH . For the special case of a monodisperse 
system, we simply set X2 = 0. 

We may also calculate the nematic order parameters, Si^m, which define the local order of the m-th segment 



Si.m ^ J dujfi,miu})P2{u} -n), (18) 

where P2('^ • n) is the second Legendre polynomial, and n is the nematic director. We define the nematic order of a 
rod as the average nematic order along the chain 



M, 

S^-JtT. S..n. (19) 



m— 1 



Effective length 

The calculation of the effective length goes as follows. We define the mean square effective length ^ as 

Mi Mi 

m— 1 m' — l 

Mi Mi 

m— 1 m' — l 

where we are summing the squares of the average length projections of all chain segments m! along the director of all 
segments m. Here, fi^m.m'i'-^,'-^') is the pair orientational distribution function (PDF) defined by 

/i,m,m'(w,a;') = / fi{n)5{uj,n ~ Lu)5{uj,n' ~ Lu')dn, (21) 



where we are integrating out all other degrees of freedom from fi{fl) except those of segments m and to'. Note that 
fi,m,m'i'^, is the probability that a chain of species i is in a configuration with the TO-th and TO'-th segment having 
orientations uj and w', simultaneously. Inserting Eq. ([ro]) into Eq. ([2T|) . and using Eqs. (fT4|) and (fT6|) . we find that 

/j,m,m'('^,w') = -^gi,m(w)exp[-i;i(w)](3i^„,™'(a;,a;')exp[-u,;(w')]gi,M-m'+i(w'), (22) 

with the same notation as before. Here Qi^m.m'i'-^T'-^') is the partial chain partition function that takes into account 
the effect of the chain segments that link segment m and segment m'. For neighboring segments m' = to + 1 we have 
Qi.m,m'(w, uj') = exjp[—Ui{uj , uj')], and for to' = to + 2 . . . M , it follows the recursion relation 

Qz,7n,ni'{i^,oj') = / da;"(3i,m,m'_i(w,w")exp[-Wi(w")]exp[-Uj(w",w')]. (23) 



By construction, each pair orientation distribution function also obeys the normalization condition 
J fi,m,m'{uj,u!')dujduj' = 1. As wc already know the ODFs, and hence qi.mi^j) and Vi{uj), from our phase diagram 
calculations, the calculation of the PDFs and Le,i is relatively straight forward. We use Le,i to calculate the diameter 
Dg^i required for rigid rods to have the same excluded volume as our flexible rods, at the same state point, obtaining 
the effective shape of the rods. 
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fd virus parameters and values 



In order to accurately describe semi-flexible rods, we must ensure that we are in the continuum limit. That is, we 
must use a sufficient number of chain segments to ensure that our results capture the physics of a continuous chain. 
We do this by checking the convergence of our results with increasing Mi at fixed Pi, by adapting k — Li /Mi. We 
shall now examine the one-component system fd virus system, and hence we drop the subscript i. In Fig.|4l we show 
the isotropic (I) and nematic (N) equation of state of a one-component fd virus system, for various M values, starting 
from the rigid-rod limit M = 1. Clearly, the isotropic branch is independent of M , while the nematic branch strongly 
depends on M. The I-N coexistence, which is represented by the jump, shows a phase transition that shifts to higher 
p and n upon increasing M, reaching a well-defined continuum limit in the pressure regime of interest for M > 15. 
Coexistence is found at HB = 29.54. 

Varying the stiffness of semi-flexible rods can have a large effect on the properties of the nematic phase at coexistence. 
Fig. [5] shows the equation of state for bio-engineered fd virus particles of various persistence lengths. We see that 
coexistence is found at much lower densities (and osmotic pressures) for stiff rods than for flexible ones, in agreement 
with Refs. 

We may also calculate the nematic order parameter of each chain segment in a rod from Eq. (fT8|) . Fig. [6] shows 
the nematic order parameter S' at a distance r € [0, L] along the rod, at isotropic-nematic coexistence for fd virus 
particles of various persistence lengths. We see that rigid rods are more ordered at coexistence than flexible rods, 
despite coexistence being found at a lower density. The chain segments in the middle of the rods are also found to be 
more ordered than the end segments, in agreement with earlier bifurcation findings in Ref. [l5| . 




pB 



FIG. 4. Convergence of equation of state in pB — IIS representation, for increasing number of chain segments M for fd virus 
parameters L — 0.88/im, D — 6.6nm and P — 2.2/^m. We find that Al — 15 is within the continuum limit for the fd virus. 
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FIG. 6. Nematic order parameter S of fd virus particles of various persistence lengths, at a distance r G [0, L] along the rod. 
The densities of the nematic phase at coexistence are: pB — 4.19, pB — 4.27, pB — 4.42, pB = 5.45, pB — 7.17 and pB = 11.30 
for P oo (rigid-rod), P — 12/im, P = 6.5^m, P = 2.2p.m, P = 1.1pm and P — 0.55/xm, respectively. 



